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First post-Newtonian (1PN) hydrostatic equations for an irrotational fluid which have 
been recently derived are solved for an incompressible star. The 1PN configurations are 
expressed as a deformation of the Newtonian irrotational Riemann ellipsoid using Lagrangian 
displacement vectors introduced by Chandrasekhar. For the 1PN solutions, we also calculate 
the luminosity of gravitational waves in the 1PN approximation using the Blanchet-Damour 
formalism. It is found that the solutions of the 1PN equations exhibit singularities at points 
where the axial ratios of semi-axes are 1 : 0.5244 : 0.6579 and 1 : 0.2374 : 0.2963, and the 
singularities seem to show that at the points, the irrotational Riemann ellipsoid is unstable 
^ ' to the deformation induced by the effect of general relativity. For stable cases [0,2/0,1 > 

ON ' 0.5244, where ai and 0,2 are the semi-major and minor axes, respectively) we find that when 

C*") . increasing the 1PN correction, the angular velocity and total angular momentum increase, 

while the total energy and luminosity of gravitational waves decrease. These 1PN solutions 
will be useful when examining the accuracy of numerical code for obtaining relativistic 
\ irrotational stars. 

We also investigate the validity of an ellipsoidal approximation, in which a 1PN solution 
is obtained assuming an ellipsoidal figure and neglecting the deformation. It is found that 
for 02/ai > 0.7, the ellipsoidal approximation gives a fairly accurate result for the energy, 
angular momentum, and angular velocity, although in the approximation we cannot find the 
I . singularities. 
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§1. Introduction 



Preparation of reliable theoretical models on the late inspiraling stage of binary 
neutron stars is one of the most important issues for gravitational wave astronomy. 
This is becauseJhey represent one promising source for gravitational wave detectors 
such as LIGO,B VIRGO,! GEO600& and TAMA.i From their signals, we will 
obtain a wide variety of physical information on neutron stars such as their mass 
and spin, if we have a theoretical template of th.em.fP In particular, a signal from 
the very late inspiraling stage just prior to merging may contain physically impor- 
tant information on neutron stars such as their radius, which will be utilized for 
determining the equation of state of neutron stars. H 

Binary neutron stars evolve due to the radiation reaction of gravitational waves, 
so that they never settle down to equilibrium states. However, the emission time 
scale will always be longer than the orbital period outside their innermost stable 
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circular orbit (ISCO), so that we may consider them to be in quasiequilibrium states 
in their inspiraling phase even near the ISCO. 

Until now, all the reliable relativistic works devoted to obtaining a quasiequi- 
librium state have been undertaken assuming a corotational velocity field, W'U 1 since 
there was no formalism to compute non-corotational solutions. As pointed out pre- 
viously, a>' however, corotation is not an adequate assumption for the velocity 
field of realistic binary neutron stars, because the effect of viscosity is negligible for 
the evolution of neutron stars in a binary system and, as a result, their velocity 
fields are expected to be irrotational (or nearly irrotational). Formalism for the 
computation of realistic quasiequilibrium states of coalescing binary neutron stars 
just prior to mexging has recently been developed by several authors in general 
relativity. t!P'E3)> t3' ES >i3 In this formalism, we have only to solve two hydrostatic 
equations as for the fluid equations. One of them is the integrated form of the Euler 
equation and the other is the Poisson equation for the velocity potential. Thus, the 
formalism seems to be very tractable for computing equilibrium configurations of 
relativistic irrotational bodies. 

In this paper, we apply the formalism for solving an incompressible, irrotational 
single star in the first post-Newtonian (1PN) approximation as a first step. Incom- 
pressible rotating stars were studied extensively at the 1PN order by Chandrasekhar 
about 30 years ago. He derived many equilibrium configurations of the Maclaurin, 
Jacobi, Dedekind ellipsoids, developing original methods (see also Ref. 

20)), and he found many interesting features of rotating stars deformed by relativistic 
effects. 

The purpose of this paper is twofold. One is to develop a formalism to obtain an 
equilibrium state of an irrotational, incompressible fluid in the 1PN approximation. 
As mentioned above, Chandrasekhar studied 1PN equilibrium stars extensively, but 
he did not do so for irrotational stars. We develop the method in this paper which 
can be applied even for irrotational binary systems. The other comes from a demand 
in performing numerical computation: Our final aim is to study general relativistic 
irrotational binary neutron stars with compressible equations of state. Thus, nu- 
merical computation is necessary. Although some promising numerical methods are 
proposed, the computation still does not seem easy. Hence, when we obtain a 

result by numerical computation, it is necessary to check its validity by comparing 
with an analytic solution. Thus, preparation of exact solutions is required. 

This paper is organized as follows. In §2, we describe the basic equations to cal- 
culate equilibrium configurations of the 1PN irrotational star. The deformation of 
the star from an ellipsoid and the angular velocity in the 1PN approximation are cal- 
culated in §3. In §4, we give the boundary conditions to determine the velocity field 
and the deformation in the 1PN approximation, and the total energy and angular 
momentum are obtained in §5. We present formalism for computation of gravita- 
tional radiation from a rotating star with arbitrary internal motion in §6. In §7, we 
present the numerical results obtained by solving the equations derived in previous 
sections. We also compare the results with those in the ellipsoidal approximation. EHP 
Section 8 is devoted to summary and discussion. 

Throughout this paper, c denotes the light velocity and we use the units in which 
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G = 1. Latin indices k, ■ ■ ■ take values 1 to 3, and 5ij denotes the Kronecker delta. 



§2. Formulation 



Non-axisymmetric equilibrium configurations with non-uniform velocity fields 
are obtained by solving the Euler, continuity, and Poisson equations, consistently. 
Since we consider an incompressible fluid, all the calculations are carried out ana- 
lytically even in the 1PN case. The procedure is as follows. 

(1) We calculate a Newtonian equilibrium state, i.e., the irrotational Riemann ellip- 
soid,!^ as a non-perturbed state. For simplicity, we consider only the case in which 
the directions of the vorticity vector and the angular velocity vector lie along X3-axis. 

(2) 1PN corrections for the velocity potential and gravitational potentials are ob- 
tained from the 1PN Poisson equations for them. 

(3) We calculate the deformation from the Newtonian ellipsoid induced by 1PN grav- 
ity using Lagrangian displacement vectors introduced by Chandrasekhar.EU Then, 
corrections of the Newtonian quantities due to the deformation of the star are esti- 
mated. 

(4) We substitute all the 1PN corrections obtained in (2) and (3) into the 1PN Euler 
and continuity equations. Then, coefficients of the Lagrangian displacement vectors 
and 1PN velocity potential are determined from the boundary conditions on the 
stellar surface. 

In this section, we calculate the Newtonian and 1PN terms which we need in the 
above procedures. In the following, we assume that the center of mass of the star is 
located at the origin of the coordinate system, and the direction of the coordinate 
axes are parallel to the principal axes of the star, whose lengths are defined as oi, «2 
and CL3. 

2.1. Hydrostatic and Poisson equations for 1PN irrotational stars 

For an irrotational fluid, therelativistic Euler equation can be integrated, and 
in the 1PN case, it is written asi3 

const = --U + \ ]T(<9 fc N ) 2 - Y, &dk<fo 
p Z k k 



1 

+- 

& 



~U + \u 2 + X - \ {t + 3U ) $>^n) 2 _ ± ( £(^n)' 



(2-1) 



where we assume that fluid is incompressible, i.e., p =const. In Eq. (2T), P, p, £ k , 
0Nj 0pn, U, X, and /3 k , respectively, denote the pressure, the density, the velocity 
field of the figure rotation (spatial component of the Killing vector), the Newtonian 
and 1PN velocity potentials, and the last three terms are the Newtonian and 1PN 
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potentials, which are derived by solving Poisson equations as 

AU = -47rp, 



AX = 4irp 



op 

2U + — + 2J2(9M 



AP k = -4npd k (j) N , 
A X = 4irpJ2(dk(l)N)x k . 



(2-2) 
(2-3) 

(2-4) 
(2-5) 



Here Pk is expressed as 



Pk 



(2-6) 



We note that using the gravitational potentials, the spacetime line element to 1PN 
order can be written as 



ds' 



-a: 



- 2 c 2 dt 2 + — Pidx { dt + 1 + — )J2 dx l dx 



2U 



where 



U 1 /U 



a = 1 + 

c 2 cf 



+ X)+0(c- 6 ). 



(2-7) 



A characteristic feature of the irrotational fluid is that the continuity equation 
reduces to a Poisson type equation for a velocity potential </),t3)'EJ) and in the 1PN 
incompressible case, it is 

X>diC; = 0, (2-9) 



where 



i (\ E( 9 ^n) 2 + 3U) - -di<h ~ Pi + ds 



• (2-10) 



In the following subsections, we obtain terms appearing in Eq. (2-1) by separately 
solving equations for each. 

2.2. Newtonian terms 

The solution of the Newtonian potential for the ellipsoidal star is 

U = 7rp(A -J2Ak4), (2-H) 
k 

whereAj... are index symbols introduced by Chandrasekhar,© and Aq is calculated 
from 



Aq = aia 2 a 3 / 
Jo 



da 



a l a2Ct3 



^{a\ + u){al + u){al + u) 
dt 

l + t)(a 2 2 + t)(al+t) 



(2-12) 
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where 02 = 02/01 and «3 = 03/01. Also, the pressure at Newtonian order is written 
as 

„2 ■ 



\ u k ' 



k "fc 

where -Po = ^P 2 ^^ denotes the pressure at the center of the star. 

The Newtonian velocity field in the inertial frame v% = ft^N is written as 

vi = I 1 + u\ = - ( 9 1 ^' fl 9 + l) (2x2, 
Vof + o^ / 

t>2 = i 2 + ^2 = ( o 2 ^ R 9 + 1^) 

Vaf + o^ J 

v 3 = 0, (2-14) 

where Ui is the velocity field in the corotating frame and (2 denotes the angular 
velocity of the figure, t and U{ are given by 

t = (-f2x 2 , (2xi, 0), (2-15) 

f a l a 02 

Ui = — Ax 2 , Axi, 

\o 2 Ol J 

■4^^, 0), (2-16) 

ti^ HQ (*1 "T (^2 / 

where yl is the angular velocity of the internal motion, fn is defined as C,/(2, where 
( = (rotu)3 denotes the vorticity in the corotating frame. For the irrotational 
ellipsoid, /r becomes —2. Thus, 0n = F a (2x\X2 and 

v i = (F a (2x2,F a f2x 1 ,0), (2-17) 

where 

F ^ a M- (2 - i8) 

2.3. 1PN terms 
2.3.1. X 

For the convenience of the calculation, we separate X into two parts as 

X = X + X V , (2-19) 
where Xq and X„ are derived from Poisson-like equations as 

AX = 4irp (2U + ^pj , (2-20) 
AY„ = 87T P J2{dk<f>N) 2 . (2-21) 
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Then, the equation for Xq becomes 



AX = Airp 
and the solution is 



X = -a U + ^2 VkDkk, 
k 



T 2 

x k . 



where 



3P 

a = 2irpA H , 

P 

7] k = 2TrpA k H 2 » 



2 \ 2 



4 (^fcfc — X! ^ fcfc 
m 

BkmXm + Bk 



2 2 



(2-22) 
(2-23) 

(2-24) 
(2-25) 

(2-26) 



and Bij... are index symbols introduced by Chandrasekhar. S Dy... is calculated 
from a Poisson equation as 



ADij... = —AirpxiXj 



The equation for X v is 



zAX„ = 8irpF 2 f2 2 (x( + 



and hence, the solution is 



X v = -2FlQ 2 (D ll +D 22 ). 



(2-27) 
(2-28) 
(2-29) 



2.3.2. k 

Substituting Vi of Newtonian order into the equations for Pi and Xi we immedi- 
ately find the solutions as 



P 1 = F a QD 2 , 
P 2 = F a flD 1 , 

^3 = 0, 

X = -2F a nD 12 , 



where 



A = irpajfAi-^A 
k 

D 12 = Trpalal(^Ai 2 - ^ A 12k x 2 ^jxix 2 . 



(2-30) 
(2-31) 
(2-32) 
(2-33) 

(2-34) 
(2-35) 
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Using the above solutions of Pk and Xi we obtain (5^- 

■h. Ti/iFJi \^a\A x - -a\A 2 - a\a\A 12 \x 2 

+ \--a\Axi + ^alA 12 + 3a\a\A n2 )->i->'2 



1 
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+ ( -~ajA 12 + -a z 2 A 22 + a\a l 2 A X22 )x 6 2 

+ ( -^' l - 1 !-'. + 7^2^23 + "T > "2-1|2:! ) •^■'"a 



(2-36) 



-a 2 A 2 — -afAi — a\a 2 A\ 2 \x\ 



+ | ~a%A 12 + ^i^n + a\a\A xvl ) y\ 
+ ( ~a^A 22 + -a\A 12 + 3ajalA 122 )x t xi 



+ ( -^2^23 + ^a?Ai 3 + af a|A 123 j »' i .r : , 



/3 3 = 7rpF a J?f2aia!^i23 - a|^23 - a^i 3 )xix 2 2; 3 . 
Then, the divergence of ^ is 

X! 5 fcAs = 6npF a f2(aj + 03)^12x1x2, 



(2-37) 
(2-38) 



(2-39) 



where we have used some relations among index symbols. 
2.3.3. (/>PN 

An explicit form of the Poisson-type equation for 0pn can be written as 



6vrpf?(Ai - A 2 ) - 2F a Q ( + — 

a 1 a 2 J p 



X\X 2 



-2Q 



(ai-a 2 2 ) Po 

1 1 



■x\x 2 . 



(2-40) 



The solution of this equation can be written up to biquadratic terms in x% as0 

0pn = (p + Q%i + rx 2 + sx^)x\x 2 + const, (2'41) 
where g, r and s satisfy the condition 



3q + 3r + s = — Q- 

a{a 2 P 



(2-42) 



*' Although we can add higher order terms which satisfy Z\0pn = 0, we negelct them for 
simplicity because we consider only biquadratic deformtion of ellipsoids in this paper. 
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2.4. Collection 

Substituting the terms derived above into Eq. (2-1), we obtain 



P 
P 



g^- n -U-5U + ^F 2 a (Ql + -^5Q 2 ){xl + xl) 



1 



V J Km 7 



const, 



(2-43) 



where J7r denotes the angular velocity of an irrotational Riemann ellipsoid and the 
jij are expressed as 



70 = T^P 



71 = 



-F 2 Q 2 R (a 2 B 1 + a 2 B 2 ) 



P P 
4^i— + ■KpA Q A 1 + A — t 
P paf 



(2-44) 



2a 2 i^R 



72 = VT/9 



+F a 2 ^(3a^! - a\A 2 - 2a 2 3 A 3 - 3ajA u ) - - -p 

a| + a 2 vt/9 

fn -R / 3P \ ( 1 \ 

4A 2 — + irpA A 2 + A -^ + 2vrpa^ 2 + — 0^22 - -£ 22 
P pa 2 V P / V 2 / 



(2-45) 



+F a 2 ^|(3a^ 2 - a 2 Ax - 2a 2 3 A 3 - 3af A 22 ) + — 

af + vt/9 



(2-46) 



73 = ^p 



4A 3 ^ + ir P A A 3 + + ^pa 2 3 A 3 + ^) (a§4 3 3 - ^33) 

-^wpalAi + ^)si 3 - ^27rpa£A 2 + B 23 



+F 2 [2 2 R (a 2 1 B 13 + a^ 23 ) 



(2-47) 



711 = t*P 



B 



113 



712 = Trp 



(2-48) 

#112 ) 



713 = 



+i (2ixpa\A 2 + ^)#m + \ (27rpa 2 3 A 3 + ^ 
F 2 Q 2 

+^-/( 3 «?^i + a i A z - 8ofAn + 5a?A m - a\a\A 

F^Q R 2a| Q R 1 
87rp + a| vt/9 

"7^ + f ) + ^ 2 + + v) H* 1 ' + s , 

+ (2irpa 2 A 2 + ^-^ (-a^m + ^#122) + ^2vrpa|A 3 + ^)#123 

FlQ\( A A alial + al) A 
2 V " afai 

-4(ai + ai)Ai2 + 3a?(a? - a^n 2 - 3a 2 (a? - al)Ai 22 ^ 

-Sfa + ^f^-^fkJ (2.49) 
4-rrp af + a 2 vt/9 af + a 2 vrp 

-7(1 + + TMlA3 + l 2 *^ 1 + T 1 ) H- 4 ' 13 + ^ 113 ) 

+\(tocp$Ai + ^)Si23 + (2^3 + ^) (-a§A 133 + ^33) 

(2-50) 



^(4A 3 - a^ 23 - 9a?A 13 + 6afA 113 ) - -^L^ 
2 af + a 2 vt/9 



722 = T^p 



F 2 Q 2 

±_a_R(± As _ a 2^ 28 _ g a 2 Ais + 6a 4 Alls ) 

af + a 2 7rp 

-A 2 ^\ + \vpA\ + (2irpc%A 2 + — ) (~a 2 2 A 222 + \b 
pa 2 2 V P / V 1 



4 

#223 



222 



+- ^2vrpa^i + -y)Bi22 + j \2iipa 2 A 3 + — 
F 2 Q 2 

+-±^(3a%A 2 + a 2 A 3 - Sa\A 22 - a\a\A 122 + ba^A 
F^ R | 2a| ^ 

87T/9 + 02 VT/9 



'2^222 J 



(2-51) 











723 = 






+ ^) 




P 


va 3 


« 2 / 



+ ^2vrpal^ 3 + ^) (-a 3 ^233 + ^#233) 

+^i(4^ 3 - a 2 A 13 - 9a 2 A 23 + 6a^ 223 ) + 

2 a( + a 2 irp 



#123 



(2-52) 
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733 = T^p 



A 3^2+ l^pAl + \2vpa\A z + 



P a 3 

+\(2itpa\Ai + 



3^ 
P 



-ai^333 + 7-B333 



P 



3-P 

2npa\A2 H 



-B233 



^ (011-0133 + 0-2^233) 



(2-53) 



<5C7 denotes the gravitational potential induced by the deformation of the figure. The 
explicit form of 5U is given in the §||. 

§3. Second tensor virial equation 

The angular velocity is determined from the second tensor virial equation,© 
which is written as 

r g 

J d 3 xp Xj —g = 0. (3-1) 

Substituting Eq. (2-43) into the above equation, we obtain 

= -<%(/! + ,577) - Wij - 5W l3 + F a {F a - 2){Q% + <W2 a )$ii( J V + 5 hj) 
+F a (F a + 2){Q 2 R + 6{2 2 )5 2i (I 2j + 5I 2j ) 
1 



+- 



2jihj + 4:~fxilmj6ii + 27i 2 (Ji22j^ii + hujfai) 

+27l3(7l33j^l + hl3j$3i) + ^I22h22j$2i 
+2^f23{h33jhi + h23jhi) + 4733^333j^3i , 



where 



n 



Wi 



■1 j 



L 'j 



d 3 xP, 

dU 

(ftxpxi- — = -2-rrpAiIij, 
dxj 



M 

T 

M 
35 



a 2 5ij, 



d 3 xpx 2 x 2 = —a 2 a 2 (l + 25i 3) , 



sn = 6 



J d 3 xP, 



5Wi. 



6 / d 3 xpx. 



dU 



5 / d 3 xpxiXj. 



M denotes the Newtonian mass: 



M 



d 3 xp. 



(3-2) 

(3-3) 
(3-4) 

(3-5) 
(3-6) 
(3-7) 
(3-8) 
(3-9) 

(3-10) 
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(577, 5Wij, and 81{j must be taken into account because the ellipsoidal figure deforms 
due to the 1PN effect. The explicit forms are given later. 

3.1. Second virial equations at Newtonian order 

First, we derive equations for the angular velocity and axial ratios of the ellipsoid 
at Newtonian order. The second tensor virial equation at Newtonian order is written 

= -<%/7 - Wij + F a {F a - 2)Q 2 R 5 li I lj + F a (F a + 2)Q 2 R 8 2i I 2j . (3-11) 



The non- vanishing components are 



n 
n 
n 



-Wu+F a (F a -2)n 2 R I u , 
-W 22 + F a (F a + 2)n R I 22 , 

-w 33 . 



(3-12) 
(3-13) 
(3-14) 



Then, we obtain 



R 



2i 



1 + "2^2 = 2 ^ 5 12, 

\a{ + a l 2 ) J 

2a\a\ 2 . 2 2 a 2 a \ 
■ 2 i 2 n R = vp(a x a 2 A X 2 - a 3 A 3 ). 



(3-15) 
(3-16) 



From these equations, we can determine the axial ratios a 2 ja\ and a 3 /a\, and the 
angular velocity. 

3.2. Deformation of the figure 

The deformation of an incompressible ellipsoid due to the 1PN gravity is written 

as 



by using the Lagrangian displacement vectors Q of J2k ®k£, 

£ fc = ^2 Sijfil^' 

ij 

where 



(ij) 

k 



t(n) 

t(12) 
t(31) 

t(32) 

t(33) 



(xi, 0, -x 3 ), 

(0, X 2 , -X 3 ), 



x^, —x\x 2 , 



0) 3^2' 



-x|x 3 



-X3X1, 0, — x| 



1 

3' 



(3-17) 

(3-18) 
(3-19) 

(3-20) 
(3-21) 
(3-22) 



Here, we consider only up to the cubic deformation and also assume triplane sym- 
metries, so that the Lagrangian displacement vectors of higher order functions in X{ 
and of even functions can be neglected. 
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Using the Lagrangian displacement vectors, we can calculate Eqs. (3-7) ~ (3-9). 
We find £3) 



(577 = 0, 



5Iij= J d^xpJ^ti-j^r&iXj), 



Shi 
SI33 



11 



5W: 



22 



5W : - 



33 



^2 (Slllll + g53i/im — S33I1133 



2 [ S\ 2 l22 — 531^1122 + ^'S'3272222 



1 



-S11I33 — S12/33 — £32/2233 + 2'S'33/3333 



SWij = Trp(-2BijSIij + a-Sy ^ Au5I u j , 



2vrp 



5n{(o?An - 2B n )I u - 0^3/33} 

+S 12 (a{A 12 I 22 - 0^13/33) 
+53i|i(o?A u - 2B n )I un - a\A l2 I ll22 



2vrp 



+S32 ^-0^1272222 — 01^13/2233 
+S 33 ^-(a 2 1 A 11 - 2B n )I 1133 + ^0^13/3333} 



5n(a|Ai 2 /n - 0^23/33) 
+S 12 \{a 2 2 A 22 - 2B 22 )I 22 ~ a\A 23 l 

+£31 j -02^12/1111 — (0^22 — 2B 22 )I\\ 22 

+5 , 32 1^ ~(alA 22 — 2B 22 )I 2222 — 0^23/2233 

+5"33 ( — O2A12/1133 + -0^23/3333 



'33 



2vrp 



Sii|o|Ai3/n - (a|A 33 - 2B 33 )/33} 

+5i 2 |o|A 2 3/22 - (a^ 33 - 2B 33 )I 33 



+S 3 i(^-a 3 Ai 3 Iim — 0^23/1122 

+5*32 1 -03^23/2222 — (0IA33 — 27?33)/2233 



(3-23) 
(3-24) 

(3-25) 

(3-26) 

(3-27) 
(3-28) 



(3-29) 



(3-30) 
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+533 j -03^13-^1133 + 2(03^-33 — 2^33)73333 



(3-31) 



3.3. Second virial equations in the 1PN approximation 

Substituting the 1PN terms into Eq. (3-2), we have the second virial equations 
in the 1PN approximation. Their components are given as follows: 
1. i = j = 1: 

12 2 2 

= 2 7l a? + y7n<4 + -ixia\al + -^a\a\ + F a (F a - 2)a\5Q 2 



+2F a (F a - 2)Ql{s lia \ + ^S 31 a\ - ^ 33 a?a|) 



-27rp Sni^An - 2a 1 B u - a^A^) + S^a^A^ - a^A^) 

+ ^f(a e 1 A n - 2afB n - a\a 2 2 A 12 ) + ^{a\a A 2 A l2 - a?^a§A 13 ) 



+^{-a\alA u + 2a\alB lx + a\a A 3 A 13 ) 



(3-32) 



2. i 



2 1 ^ 2 2 , 12 4 2 



= 2-f 2 a 2 + -7i 2 aia 2 + y722a 2 + -723a 2 a 3 + i^C^a + 2)a 2 5Q 
+2F a (F a + 2)r?^5 12 a| - ^S 31 aja 2 + ^S 32 a^ 



-2tt P 



S n (a\alA 12 - a 2 a 3 A 23 ) + S 12 (a 2 A 22 - 2a 2 B 22 - a 2 a A 3 A 



,2„2 



2„2 



23 I 



(3-33) 



+^(0^22 - 2a|5 2 2 - aja^A^) + ^y-i-ajalalAn + ala 3 A 23 ) 
S 

+—{a\alAi 2 - a\a\A 22 + 2a\a 2 2 B 22 ) 

3. i = j = 3; 

= 273^ + ^7l30?a| + ^723 02«3 + y 733^ 



-27T/9 



Sn(afa§Ai 3 - a|A 33 + 2a§ J B 33 ) + S 12 {aiai,A 23 - a%A 33 + 2a z 3 B 33 ) 
+—(03^33 - 203/^33 - 0103^13) + — (a\a\A 13 - a\ala 2 3 A 23 ) 



S 

+ ^r( a 2 a 3^23 - a2«3^33 + 2a|ci3.B 33 ) 



(3-34) 



Subtracting the (2, 2) component of the 1PN second virial equation from the 
(1, 1) component, we have the equation for Sf2 2 as 

= 271a? - 272^ + -(671101 + 7i3a?a| - 6722^ - 723ai«3) 



14 



-{a\ - a\) 



1 + 



2a\(i2 
a\ + a 2 



-27TP 



SllTl + S , 12T2 + ^-5 < 3i(ti — T2) + -yS l 32T2 ^S^Tl 



(3-35) 



where we use the equilibrium equations of Newtonian order and define new variables 

as 

A A „2 2 a _ 2 2 A , „2„2 A __ 2 



n = 30^11 - 0^2^12 - a^A^ + a 2 a 3 A 2 3 - 2a 3 A 3 , 
r 2 = a\a\Ai2 - 3a 2 A 22 - a? 0^13 + a\a 2 A 23 + 203^3. 



Then, Sf2 is expressed as 



5Q 2 = 



irp(a{ - a%)B 12 



(3-36) 
(3-37) 

71 a\ - 72CI2 + ^(671101 + 713a? a| - 672203 - 723a|a|) 



^2 ^2 

-7rp( 5nTi + S 12 t 2 + -y<S3i(ri - r 2 ) + y5 3 2r 2 - y^ri 



(3-38) 



By calculating (1, 1) — (3, 3) and (2, 2) — (3, 3), we can find other expressions of Sf2 2 
as 

= 2 7l a? - 273a! + ^(frynaf + Juajaj - 723^ - 6 733 a|) + F a (F a - 2)a\5Q 2 



-27rp 

2 



Qry (X<2 

Sn(n - r 3 ) - Si 2 t 3 + —SziTi - —S 32 t 3 - —S 33 (t\ - r 3 ) 



-2np 



= 272^ - 27303 + ^(7120^2 + 672202 - 713a? a| - 673303) + F a {F a + 2)o2^J7 

2 



, (3-39) 

2rn2 



-5*11X3 - S\ 2 (t 2 + T 3 ) + -yS 3 \T 2 - -yS 32 (r 2 + T 3 ) + -yS 33 T 3 



, (340) 



where 



r 3 = a\a 2 3 A l3 - 303^33 - a\a 2 2 A l2 + a 2 2 a 2 3 A 23 + 2a|^ 3 - (3-41) 



Eliminating 8fi 2 from Eq. (3-39) by using Eq. (3-40), we obtain the identity 

1 



= (3a? + a 2 2 )a\ 
+{a\ + 3al)a\ 
— irp(3a 2 + a 2 )a 2 

—irp(a 2 + 3a 2 )a 2 



7i«i - 730 3 + -(6711a! + 7120^2 - 723a 2 a 3 ~ 673303) 



7202 - 73a 3 + ^(71201^2 + 6722a 2 - 713^103 - 673303) 



Su(ti - T 3 ) - Si 2 T 3 + -yS 3 lTi 



-S32T3 ^-S 33 (ti - T 3 ) 



-Sl\T 3 - Si 2 (T 2 + T 3 ) + — S 3 \T 2 

a\ a 2 
— yS 32 (T 2 + r 3 ) + ^533X3 



(3-42) 



Equation (3-42) is one of the equations to determine p, q, r, s, and S. 
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§4. Boundary conditions 

In this section we derive equations to determine the coefficients of the 1PN veloc- 
ity potential and the deformation of the figure, p, q, r, s, and SV,-, from the boundary 
conditions on the stellar surface. The stellar surfaces of the Riemann ellipsoid and 
its deformed figure are expressed as Sr(x) = and S(x) = 0, respectively, where 

„2 



Sr(x) 
S(x) 



(4-1) 



E 

i 



XT 



dS 



R 



J dxj 



Sr(x) 







Sii [ -4- 


4\ 






~ 4) 



+ s 



12 



44 



™2 
J. 2 



+ S; 



311 3a? 



2 9 

x 1 x 2 









44\ 


V 3al 


a\ )\ 



(4-2) 



The boundary conditions for the continuity equation and integrated Euler equation 
are, respectively, 



dS 



( A ) E C ^ = ° on S(x)=0, 



j 



on S(x) = 0. 



(4-3) 
(4-4) 



The condition (A) comes from the constraint that the normal component of the 
velocity on the surface vanishes. The condition (B) determines the stellar surface. 
Eqs. (4-3) and (4-4) must hold to order of c~ 2 . 



4.1. Condition (A) 

Substituting Cj given by Eq. ( |2-10| ) into Eq. (4-3), we obtain an equation of 
1PN order as 

= xix 2 {cq + c\x\ + c 2 x\ + C3X3), (4-5) 

where Ck (k = 0,1,2,3) are functions of a^, p, q, r, s, and Sij. Equation Q4-5|) must 
be satisfied on Sr. Then, we obtain the following three equations for p, q, r, s and 

Si 



c 



(A) 



4 2 P ' 2" 2 9 ' > 



aTa 



1"2 

3irpf2 ( — £ 
npF a Q 



1"2 

-4i 



af + a 2 



1 



+ 5; 



12 J 



:5 ; 



31 



1 



(A) 



+ 2a? 



.a 



1 )B n + 

4f2 



1 

2' 



3 U 1 



a\a\ 



3af + a^ 

2„4 P ^ 1272 r T2 



«1«2 



aia 2 



+ o 2 L«2 



7-f -1 B 12 + 4(^ + a^)B 112 
a i J 

-(— 5*11 + S12) H — \Ssi + - 



0,(4-6) 



S: 



°2 



-■',2 
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+ 



37T/9i7 

npF a Q 



\ a% J af 



2 4 



2„2 
a 2 



-F 2 f2 3 — 

2 a ~2„2 



a? 2 



2a| 



a 



(A) _ aj + a\ 



P + 



1<* 

1 

«2 

1 



1)B 12 + 7 



i 



2 2 ' 2 
J l a 2 a 3 



S + 



4^ 



I )B 2 2 + Ka\+ai)B 122 
1 



0,(4-7) 



+ a 2 



4 



(S u + S 



12 



-|S32 + 5*33 
«3 



-3irpf2[ — A3 



n 



2„2 



ara: 



l"-2 



2a§ 



7 «1 

2 

«2 



1)513+17 



1 5 23 + 4(^+^)5 



123 



We can make two equations which do not contain Sn, S12 and p as 



C [ A) a\ 



C2 ^a 2 
C3 



0. 
0. 



0.(4-8) 



(4-9) 
(4-10) 



Thus, Eqs. (4-9) and (4-10) depend only on the coefficients q, r, s, S31, S32 and 
S33; i.e., the coefficients for biquadratic deformation. This is an important feature 
in determining these coefficients. 

4.2. Condition (B) 

To obtain Q consistently, we must calculate the correction of the self-gravity 
potential due to the deformation of the figure SU, which is given byta' 



(4-11) 



where 



\x — x'\ 



(4-12) 



Substituting the Lagrangian displacement vectors into Eq. (4-12), we have 



5U^ = 


-D 1A + 5 3i3 , 


(4-13) 


SU™ = 


--E>2,2 + ^3,3, 


(4-14) 


sum = 


~g5m ( i + 5u2,2, 


(4-15) 


suw = 


— -5222,2 + 5223,3 j 


(4-16) 


6UW = 


— g5333 ; 3 + 5133,1, 


(4-17) 



where for i 7^ j, 
Dm = Trp 
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3 



I m 



(4-18) 



ijlm%l X m ) Xj 



(4-19) 



Substituting 5U and the equation S = into Eq. (2-43), we obtain 



PJs 



2Pq 

pc 2 



t L 1 tL 



Sl1 1 ~¥ ~~ 7§ + 5,12 7§ ~~ dr ) + ^31 ( tt-T 



as a? 



3a 2 



3a| 



t» 2 ™ 2 
x l x 2 



5'ii(^i,i — ^3,3) + Syi{Pi,2 — B>zfi) + -531 ^-Dm,i — Du2,2^ 

+Sm(y^B>222,2 — ^223,3^ + S33 ^"^133,1 + ^-^333,3^ 



-JL( x l + x l)+F a (xt-xi) 



~~2 (j° + H W + 51 "flmxfaV) + const 

Qixj + Q^ 2 , + Q3Z3 + Qnxf + Q 2 2^2 + Q33X3 



+Q12XXX2 + Ql3^1^3 + Q23^2 X 3 + Const 



(4-20) 



where 



F (la 2 \ 

Qi = -71 - -j{F a - 2)5Q 2 - SnTrpi^Az - 3a 2 A n + a 2 3 A 13 j 

+S 12 TTp(alA 12 - 03^13) 

-S , 3i7r / 9a 2 ^aiAin - a\a\A X v2 - ^afBm + -a|B U 2^ 

a 2 a 2 
+S 3 27TPy (^2^122 - 03A23) + 5 33 7rp^(a|Si33 - 3o?Bn 3 ), 

Q2 = -72 - + 2)5Q 2 + 5n7rp(a 2 A 12 - a 2 A 23 ) 

/2a 2 \ a 2 

-Si 2 7r/9^-^^3 - 3a^ 22 + a|A 23 J + S 31 irp-±(a\B lvl - 3o|Si 22 ) 

-^TTpa 2 ^a|A 22 2 - ai a 3^223 - ^02-B 2 22 + ^03^223^ 



(4-21) 
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+3337^(^5233 - Cl^iaa), (4-22) 

Qs = -73 + S u irp(2A 3 + ajA 13 - 3a 2 3 A 33 ) + S 12 irp(2A 3 + a 2 2 A 23 - 3a 2 3 A 33 ) 

CI 2 CL 2 

+S 3 iTTp^-(a\Bii 3 - alB 123 ) + S 32 it p^{a\B 223 ~ 3a§B 233 ) 

-S 33 -npal (a 3 A 333 - a\alAi 33 - ^a\B 333 + ^a 2 B 133 ^J , (4-23) 

/2a 2 5 5 1 \ 

Qn = -7ii - S 3 i-k pa\[—^A 3 - -af^mi + ajalA ul2 + -afBim - -ajBuu J 

a 2 a 2 
-S 32 TT P ^(alB U22 - a 2 3 B 1123 ) - S 33 Trp^-(a 2 3 B 1133 - 5a 2 5 m3 ), (4-24) 

a 2 

Q22 = -722 - S 3 iir p— {a\B n 22 - ba\Bi 222 ) 

—S 32 -Kpa\ ^^-|A 3 — -a 2 A 2222 + a 2 a 2 A 2223 + -^a\B 2222 — -a 2 i?2223^ 
a 2 

-S 33 irp-2-(a 2 3 B 2233 - a\B l223 ), (4-25) 

a 2 CL 2 

Q33 = -733 - S 3 iTTp^(alB n33 - alBi 233 ) - S 32 i:p^{a\B 2233 - 5a 2 -B 2 333) 

/ 2 5 5 1 \ 

-S 33 i\pa 2 3 [^A 3 - -a|A 33 33 + afalA l333 + -a|fi 3333 - -a 2 B l333 j, 

(4-26) 

Q12 = -712 + S 3 \-Kpa\ ( -^\A 3 + a\Am 2 - 3a 2 a 2 ^n22 - ^ai-Bni2 + \a 2 2 B\\ 22 

+S 32 TTpa2 ^2^1222 — 02 a 3^1223 — -^2^1222 + -G 3 -Bl223^ 

S 33 7T P ^ (a\B 1233 - 3a 2 B 1123 ^ , (4-27) 
Q13 = -713 + S 31 Trpal (^afA 1113 - ajalA 1123 - ^ajB 1113 + ^a 2 Sn 23 ^ 

—S 32 TTp-^-^a 2 Bi 223 — 3a 2 i?i2 33 ^ 

/ 2 3 3 \ 

+S 33 TTpal\^A 3 + a|^i 333 - 3a 2 a|An 33 - -a 2 3 Bi 333 + -a\Bu 33 j , (4-28) 

Q23 = -723 - S 3 i7rp^ (alB U23 - 3alB l223 ^j 

( 2 3 3 \ 

-^A 3 + a\A 2223 - 3a\a\A 2233 - -03-82223 + -a§^ 2 233j 

+S 33 npal ^a|A 2333 - a 2 a 2 ^i 233 - ^a 2 3 B 2333 + ^ajBi 233 ^ . (4-29) 
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Since (P/p)s must vanish on the surface, we have five conditions, 

afQn + 4Q 22 - a\a 2 2 Q l2 = 0, (4-30) 

4Q 2 2 + 4Qs3 ~ a 2 2 a 2 3 Q 23 = 0, (4-31) 

afe + 4Qu - a 2 3 alQ 13 = 0, (4-32) 

<4Q u - 4Q 22 + of Qi - a\Q 2 = 0, (4-33) 

4Q33 ~ 4Qn + 4Qs ~ 4Qi = 0. (4-34) 

When we substitute 5f2 2 into Eqs. (4-33) and (4-34), we have two identical 
equations: 

(4Q33 + 4Qn - 44Q13) - (4Q22 + 4Q33 - 44Q23) 

= 7(afQn - 4Q 22 + a\Q x - a 2 2 Q 2 ), (4-35) 

{4Q22 + 4Q33 ~ 44Q23) ~ (4Qll + 4Q22 ~ 44Q12) 

= 7(4Q33 ~ 4Qn + 4Q3 ~ a\Qx). (4-36) 

This implies that Eqs. (4-33) and (4-34) can be derived from Eqs. (4-30) ~ (4-32), 
and only three of five conditions (4-30) ~ (4-34) are available. Here, we choose Eqs. 
(4-30) ~ (4-32) as the equations to be solved. Then, the three equations do not 
contain p, Su and S\ 2 , so that they depend only on q, r, s, S31, S32 and 533. This 
means that if we use Eqs. (4-30) ~ (4-32) together with Eqs. (2-42), (4-9) and (4-10), 
we can determine the six coefficients. 

The question arises when we determine the rest three coefficients p, Su and 
S12 because we have only two equations (3-42) and (4-6) (or (4-7) or (4-8)) for the 
three variables. There is one degree of freedom. This is because there is no unique 
definition of the 1PN solution as the counterpart of a Newtonian solution, as pointed 



out by ChandrasekharQ)© and Bardeen.ll (Or, we may say that we have a degree 
to fix the PN correction of the angular velocity {5Q) or axial ratio (Su or S\ 2 ) in 
obtaining a PN equilibrium solution.) There are many possible ways to compare 
a 1PN solution with a Newtonian solution, and to fix the 1PN counterpart to a 
Newtonian solution, we must impose an additional condition. (For example, Bardeen 
has proposed that we should compare two solutions fixing the angular momentum 
and the baryon number. E3) 

Since the method for connecting a 1PN solution to a Newtonian solution can be 
arbitrarily chosen, in this paper, we simply give the condition as 

S n = S 12 , (4-37) 

i.e., we fix the PN correction of the axial ratio. The reason why we choose this con- 
dition is simply that Sn and S12 in Eqs. (4-6) ~ (4-8) are erased, and manipulations 



are simplified. Then, we have three equations (3-42), (4-6), and ( 4-37 ) for solving p 



Sn, and S\ 2 . 5£2 is determined from Eq. (3.38) if we obtain these coefficients. 



§5. The total energy and angular momentum 
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5.1. Total energy 

The total energy is written as 



E = / d 3 xp 



V " J " 1 ^ 5 .4 , 5 „,2rr , 1 o i , P „,2 5 r r 2\ , ^/ -4 



C/ + -w 4 + -u 2 f/ + - Au* + -v 2 - -f/ 2 + Ofc" 4 ) 

2 2 c 2 \8 2 T p 2 ' 



1 



En + Spn + 0(tf 



where 



En 



£„4 



E 



Pv 2 /p 



Ejj2 



10 



Epn = En^pn + -E 1 ^ + ^v 2 ^ + Ep. v i + E Pv ij p + £' f/ 2, 
En^pn = -^n ^F 2 ^ 2 (a 2 , + a 2 J + 47rp. 

1 



(5-1) 

(5-2) 
(5-3) 



+ 



-j Oi* a 3* ^33 



10 



+^F a W(aL + a 2 J 



iV/^ _. _ , o 1 \ ( ^ 2 ^2 2 "\ 



2M* 
H — — Trp 



+ ^aL 5 '3i(aL^i* - a L^2*) + ^^^(aL^* - <4* A 3*) 



- 



3^o* («i* + ~ + a2*^2*} 

ai*( a L - 7a 2 u )Au + a 2 .* (a 2 * - 7a^)A 2 * 



-20^02* ( a L + o|*)^i2* 
~35 — ~p a R ^ 

-^(vrp) 2 (llA 2 , + atAl + a\A + a|^ 2 J. 



(5-4) 
(5-5) 
(5-6) 

(5-7) 
(5-8) 
(5-9) 



In the above equations, we have used Eq. (3-38) for 5fi 2 and the conserved mass 



= J d 3 xp 
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M 



1 + ^^F 2 a Q 2 R (a\ + a 2 2 ) + yvrM 



(5-10) 



instead of the Newtonian mass M. Also we substitute the mean radius of the star 
defined by the conserved mass as 



a* 



( M* \ 1/3 
Wp/3j 



(q 2 q 3 ) 1/3 



ai 



1 + 



i r i 



(5-11) 



and a;* defined by 01*02*03* = a*, <22*/ a i* = ^2/01 and ci3*/ai* = 03/01. Aijk...* 
and -Bjjfc...* are computed by using a2*7ai* and 03*/ai*. In Eq. (5-4) we have used 
some relations among index symbols.© 

5.2. Total angular momentum 

The total angular momentum is written as 



J 



d xp 



1 



,{ 1+ ^(> +w+ ^} + ! +0(c - 



</n + -t«7pn, 



where 



M, 



x 2 v\ + X1V2, 
2 



5 

JpN = </n^PN + ^ %v 2 + J VlpU + J^p/p + 



(5-12) 

(5-13) 
(5-14) 

(5-15) 
(5-16) 



- Jn 
5 



2I2r ja^Sii - a 2 2jf Si2 + ^a^5 3 i(a^ + a 



a u a 3 ^S 33 \ + 5Q{a u - a 2 *) 



2 1 
2* J 



^a 2 *5 32 



Jv v U 



J. 



VtpP/p 



J, 



7 

35 
2M*P 

35 p 
2M* 

35 



3A *(ai* - a 2 J - a u ^i* + a 2 *^4 2 * 
F a j?p(ai* - 02*), 



(5-17) 

(5-18) 
(5-19) 
(5-20) 



TrpF a f2 R (-al 1 ,a^A u - 7a^A u + a^*^* 



2 „2 



22 



+7a^A 2 * - 2a u alt,Ai2* + 20^02*^12*), ( 5-2 l) 

and 

M = |£. (5.22) 

§6. Gravitational radiation 

In this section, we derive an equation for the luminosity of the gravitational 
radiation from a rotating star. 

6.1. Frames 

In previous sections, we used the coordinate system {x^ associated with the 
principal axis of the ellipsoid, which is considered as the rotating frame. In these 
coordinates, the estimation of multipole moments is quite simple. On the other hand, 
it is convenient to use the inertial frame {Xi} for the estimation of the luminosity 
of gravitational waves. These coordinate systems are related to each other as 

X\ = x\ cos fit — X2 sin fit, 
X2 = x\ sin fit + X2 cos fit, 

X 3 = x 3 . (6-1) 
In the inertial frame, the Newtonian internal velocity is 

Vi = (-OFiXz, fiF 2 X u 0), (6-2) 
where we define Fi for i = 1, 2 as 

F i = 1 + -rr-r/fl- ( 6 ' 3 ) 

af + CL2 

Note that this expression of the Newtonian velocity is different from Eq. ( |2-17| ) 
because we do not restrict our attention to the irrotational case in this section. If 
one takes —2 for Eq. ( |6-2| ) reduces to Eq. ( [2-17 ). 

6.2. Multipole moments 

In the 1PN approximation, the energy loss due to gravitational radiation is 
written asE-T 



dE _ _J_ 
dt ~ oc^ 



(6-4) 



where (n) denotes that n time derivatives are taken, and My, and Sij are 
written as 

My = Kij + Vij + Uij + Pij + Yij + /.'; ; . (6-5) 
K i:j = [ d 3 XpX tJ , (6-6) 



23 



Yi 



2 / d'XpUXij, 

?, ! tPxpXij, 

14 eft 2 
20 d 



cPXpX'Xij, 

Rij ^ ~21(ft j d3X P^2 V k X ijk, 



Mijk = I d XpX ijk , 



fc,z 



(6-7) 
(6-8) 
(6-9) 

(6-10) 
(6-11) 

(6-12) 
(6-13) 



where e^u denotes the spatial Levi-Civita symbol, and both 'hat' and < > denote 
the symmetric traceless part; that is, 



Xij = XiXj — —5ijX 2 , 

Xijk = XiXjXk — -X 2 (5ijXk + 5jkXi + SkiXj), 
A<ij> = - (Aij + Aji) — -5ij ^2 Akk- 

k 



(6-14) 
(6-15) 
(6-16) 



First, we consider the mass multipole moments associated with the principal 
axis of the ellipsoid. We obtain mass quadrupole moments to 1PN order as 



In = / d 3 xpxf 
Jv 

Air 3 
= — pa^a?, 
15 

I 22 = / d?xpx\ 
Jv 

47T o 
= TTP a ± a 2 a 3 

15 



1 + ^2 ( S ' 11 + y Ssi ~ ~f Ss3 



+ 0(c-% (6-17) 



1 + ~9 12 — ^"--'31 + 



:-!2 



+ 0(c~ 4 ), (6-18) 



'12 



/ d?xpxiX2 = 0. 



(6-19) 



These moments are related to the reduced moments in the inertial frame as 



Ku = - {In - I22) cos 2fit + const, 



K- 



22 



-Ku + const, 



1 



^12 = 2^ 11 ~~ ^ 22 ) s i n 2^2i + const. 



(6-20) 
(6-21) 
(6-22) 
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6.3. Energy loss 

The quadrupole moments relevant to the energy loss are obtained as 



11 



P 



12 



R 



12 



V n - V 22 = (T 
Via = ^ 2 

C/ll - f^22 = 27T/9 

f/i2 = vrp 

-P22 



--(/mi + / 2222 - 6/ii2 2 )(^i 2 - F 2 ) cos Afit 
+ (/mi - I 2 222)(i ? i 2 + Fl ) cos 2ftt + const, 
--(/mi + / 22 22 - 6/1122X/? - F 2 2 ) sin 4^ 



+-(/im - / 2 222)(i ?2 + F|) sin 2fit 



Mhi -h2)-J2 Mhm - hm) 

1 

A {hi - I22) - J2 Mhm - I2211) 



3Pq 
P 

2p 
2 



(hi - I22) - X! — (hiu - I2211) 
1 a i 

(hi - I22) - X! — (hiu - I2211) 



1 a i 



+ const, 
cos 2 fit, 

sin 2 fit, 

cos 2 fit, 

sin 2 fit, 



Yu - Y 2 2 = - ^fi 2 (hm - / 2 222 + /ll33 - /2233) COS 2 fit, 
Y\2 = -y^ 2 (/llll - /2222 + /ll33 - /2233) sin2J?t, 



R11 — R22 — & l 



20 
21 



(Fi - F 2 )(/ m i + / 22 22 - 6/1122) cos 4^ 



fi 2 



~21^ Fl + F ^( Illu ~ 12222 + /l133 ~~ 12233 ) cos 2Qt 
^(Fi - F 2 )(/mi + / 222 2 - 6/1122) sin Afit 



- — (Fi + F 2 )(/nii - /2222 + /1133 - /2233) sin 2fit 



(6-23) 

(6-24) 
(6-25) 

(6-26) 

(6-27) 

(6-28) 

(6-29) 
(6-30) 

,(6-31) 
. (6-32) 



The higher multipole moments relevant to gravitational radiation in the 1PN 
approximation vanish due to the symmetry, i.e., 



M ljk =0(c- 2 ), 
Sij =0(c- 2 ). 



(6-33) 
(6-34) 



In total, we obtain 



{hi ~ I22) + ^{2(/im - / 2 22 2 )(F 1 2 + F 

+47T / 9^o(/n - / 22 ) -J2 A l(hlll ~ h2ll) 



l)fi 2 
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6-Po / 1 \ 

H ^11 — 122 — }] — % (hia — 12211) 

p V Y a ' ' 



4 16 

= + 21 (-^1 + ^2) ) (-^llll — ^2222 + A133 — ^2233)^ 



+0(c" 



(6-35) 



where temporal averaging has been done. 

Using some algebraic relations among la and Ian, the formula for the energy 
loss can be expressed in terms of (p, f2, a±, a2, Q>3: $iii ■ ■ ■ > 5*33! -^Oi /i?) 



"eft 



32 /4vr x 2 



5c 5 V 15 7 



( / oaia 2 a 3 ) 2 (af - a 2 .).!? 6 
2 



(a 2 -a 2 ) + |{2a 2 (5 11 + ^ 31 -|5 3 3 

-2a 2 (s 12 - |s 31 + |s 32 ) 
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+0i - aa)( 2 T 7r/5j4 ° + W a 2 """ a 2 JF ^ fl ' Ql ' 0,2 ' " 3 ' 



(6-36) 



where we define ^(/r, ai, 02, a 3 ) as 

ox, 02, 03) =(63/| + 102/b + ll)(af + 4) 

-(49/1 " 92/ij - 22)a\a\ - (8f R + 22)a^(a 2 + a 2 ). (6-37) 

6.4. The Energy loss rate for irrotational Riemann ellipsoids 
For a 1PN irrotational ellipsoid (/# = — 2),0 we obtain 



d£ 32 /4vr\ 
"eft ~~5c5 [teJ 



(paia 2 a 3 ) (a l - a 2 )Q 



(a 2 - al) + I [2a 2 (s n + ^5 31 - |s 33 ) - 2a 2 (s 12 - ^5 31 + |s 32 ) 
+(a x -a 2 )j — vrpA) 

59(o| + a|) - 358a? a 2 , - 6a|(a 2 + a\ 



1 

+ 147 

We can rewrite Eq. (6-38) by using the conserved mass as 



a\ + a\ 



dE 
~dt 



where 



'dE\ 



32M 2 2 



dE\ 
~dt)i 



+ 



1 /d£7 



rtt ) 



PN 



(6-38) 



(6-39) 



(6-40) 
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'dE\ 64M2 2 6 



dt J PN 125c 5 

x 



2 ( o , a l* o a 3* o \ r, 2 f o a l* o i °2* 



^1* I "11 + — "31 ^-"33 I — 2a 2 * I Oi2 ^-"31 + — "32 

2 2 x/3(5^ 2 46 1 2r >2 / 2 , 2 \ 

+ (Ol« - G 2*)j ~ ZY*P A <>* ~ Q F a n R( a U + a 2*) 

+T7= ( 59(af„ + of.) - 358a^a^ - 6a^(a^ + a|,)) 2^2 j 

14 ' \ / a l* ~r a 2* J 

(6-41) 

Note that in Eq. (6-39) there are no terms dependent on the 1PN velocity, since the 
quadrupole moment My in Eq. does not depend on the velocity. 

§7. Numerical results 

7.1. Normalization 

First of all, we change to non-dimensional parameters according to 
p ~ _ 1 ~ _ r ~ _ s 



P 



■KpAaf TTpA Tip A Tip A ' 



a _ Six ~ _ S12 a _ £31 $ — S32 d — ^33 ,„ ,s 
"11 = o' = 2' "31 = i "32 = , "33 = , U'-LJ 

7rpaf 7rpaf up np up 



y / 7rp' (7rpai) 2 ' (7r/j) 3 / 2 a 2 ^ 



and also 



Using these parameters with the relation S\\ = S12, we can rewrite equations (2-42), 
(3-42), (4-6) ~ (4-8), and (4-30) ~ (4-32) into the non-dimension forms. From these 
eight equations, we can determine eight variables, p, q, f, s, Sn, S31, S32 and 533. 

After determination of the variables, we can calculate the energy and angu- 
lar momentum of the star from Eqs. (5-1) and (5-12) and also the luminosity of 
gravitational radiation from Eq. (6-38). In the numerical calculations these values 
must be normalized. The normalized energy, angular momentum and luminosity of 
gravitational radiation are written as 

~ _ E 

E = TToi — = -En + C s £pn, (7-3) 
Mi j a* 

J = 7773 — 7777 = + C s Jpn, (7-4) 



(MX)V2 



dt \dt)/\o>al) ydt) N + Cs \dt)p N , ^'^ 



where C s is the compactness parameter defined as 



C s ee (7-6) 

<ra* 
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In the PN approximation, we assume C s <C 1. 

7.2. Ellipsoidal approximation 

The ellipsoidal approximation, in which the equilibrium configuration is assumed 
to have an ellipsoidal shape, is usefulin order to investigate features of stars or binary 
systems in the 1PN approximation, c3'E3)'E3) because the whole calculation is done 
by setting Sij = 0, and hence is very simplified. The ellipsoidal approximation 
gives an exact solution for the rotating incompressible star in the Newtonian theory. 
However, it gives only an approximate solution for the 1PN case. If the ellipsoidal 
approximation is really an excellent approximation, it becomes a robust method for 
the study of the 1PN effects. This is motivation for our investigation of the validity 
of the ellipsoidal approximation. 

Table I. The Newtonian and 1PN angular velocity, energy, angular momentum, and luminosity 
of gravitational radiation along the equilibrium sequence of the irrotational Riemann S-type 
ellipsoid. 



02/al 


az/ai 




nl 


-En 




Jn 




sn' 2 


-Epn 






(dE/dt) N 


{dE/dt) PN 


0.99 





9950 





2667 


-0 


6000 


1 


807(-5) 


8 


050(-2) 


-4.284(-2) 


2 


957(-5) 


-8 


275 


-7) 


5 


153( 


-6) 


0.95 





9741 





2667 


-0 


5998 


4 


707(-4) 


7 


837(-2) 


-4.235(-2) 


7 


312 (-4) 


-2 


160 


-5) 


1 


341 ( 


-4) 


0.90 





9464 





2669 


-0 


5991 


1 


986 (-3) 


7 


776(-2) 


-4.071(-2) 


2 


879 (-3) 


-9 


184 


-5) 


5 


660 ( 


-4) 


0.85 





9167 





2673 


-0 


5979 


4 


728 (-3) 


7 


936(-2) 


-3.775(-2) 


6 


359(-3) 


-2 


214 


-4) 


1 


348 ( 


-3) 


0.80 





8848 





2677 


-0 


5960 


8 


917(-3) 


8 


301(-2) 


-3.324(-2) 


1 


106(-2) 


-4 


258 


-4) 


2 


546 ( 


-3) 


0.75 





8505 





2683 


-0 


5934 


1 


483(-2) 


8 


841 (-2) 


-2.694(-2) 


1 


682(-2) 


-7 


269 


-4) 


4 


244 ( 


-3) 


0.70 





8136 





2688 


-0 


5898 


2 


282(-2) 


9 


515(-2) 


-1.861(-2) 


2 


338(-2) 


-1 


156 


-3) 


6 


544 ( 


-3) 


0.65 





7736 





2692 


-0 


5851 


3 


332(-2) 





1026 


-8.211(-3) 


3 


006(-2) 


-1 


757 


-3) 


9 


594 ( 


-3) 


0.60 





7304 





2691 


-0 


5789 


4 


691(-2) 





1102 


3.220(-3) 


3 


425(-2) 


-2 


591 


-3) 


1 


370 ( 


-2) 


0.55 





6835 





2683 


-0 


5711 


6 


433(-2) 





1188 


1.989(-3) 


1 


035(-2) 


-3 


739 


-3) 


2 


110( 


-2) 


f 0.5244 





6579 





2673 


-0 


5663 


7 


505(-2) 




00 


00 




00 


-4 


481 


-3) 




CO 




0.50 





6325 





2660 


-0 


5612 


8 


660(-2) 





1111 


9.301(-2) 





1499 


-5 


305 


-3) 


1 


530 ( 


-2) 


0.45 





5772 





2614 


-0 


5485 





1151 





1117 


0.1021 





1357 


-7 


403 


-3) 


2 


337( 


-2) 


0.40 





5174 





2534 


-0 


5324 





1516 





1025 


0.1379 





1699 


-1 


012 


-2) 


2 


763 ( 


-2) 


0.35 





4531 





2403 


-0 


5121 





1988 


8 


362(-2) 


0.1895 





2423 


-1 


345 


-2) 


2 


898( 


-2) 


0.30 





3850 





2206 


-0 


4864 





2606 


5 


094(-2) 


0.2762 





4161 


-1 


707 


-2) 


2 


507( 


-2) 


0.25 





3143 





1927 


-0 


4540 





3428 


-7 


291(-2) 


0.7824 


1 


632 


-2 


018 


-2) 


-3 


537( 


-3) 


t 0.2374 





2963 





1842 


-0 


4447 





3678 




00 


CO 




00 


2 


071 


-2) 




CO 




0.20 





2430 





1559 


-0 


4135 





4551 


6 


538(-2) 


5. 760 (-3) 


1 


493(-3) 


-2 


127 


-2) 


2 


552( 


-2) 


0.15 





1737 





1117 


-0 


3630 





6154 


2 


253(-2) 


0.1183 





8808 


-1 


861 


-2) 


1 


067( 


-2) 


0.10 





1094 


6 


454(-2) 


-0 


2994 





8666 


5 


964(-3) 


0.1163 


2 


973 


-1 


172 


-2) 


2 


893 ( 


-3) 


0.05 


5 


167(-2) 


2 


252(-2) 


-0 


2137 


1 


370 


5 


969 (-4) 


6.554(-2) 


14.44 


-3 


460 


-3) 


9 


916i 


-5) 



In the ellipsoidal approximation, a solution is obtained by setting S%j = in all 
the equations. After we set Sij = 0, we can calculate the velocity field in the 1PN 
approximation from Eqs. (2-42) and (4-6) ~ (4-8). We note that in the ellipsoidal 
approximation, the boundary condition (4-4) is not satisfied. 

7.3. Results 

We give the results in Tables I, II and III. In Table I, the angular velocity, energy, 
angular momentum, and luminosity of gravitational radiation of Newtonian and 1PN 
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orders are presented along the Newtonian sequence. Coefficients of the 1PN velocity 
potential and the Lagrangian displacement vectors are shown in Table II. In these 
tables, f and X denote singularities at which all the 1PN terms diverge. We believe 
that at those points, some instabilities which are concerned with the fourth order 
harmonics of the Riemann ellipsoid are induced by relativistic corrections. 

Table II. Coefficients of the 1PN velocity potential and the Lagrangian displacement vectors shown 
along the equilibrium sequence of the irrotational Riemann S-type ellipsoid. 



a 2 /a 1 


P 


q 


r 


s 


S n 


031 


532 


533 


0.99 


3.337(-2) 


-1.893(-3) 


-1.936(-3) 


-1.914(-3) 


-9.723(-6) 


4.082(-5) 


8.268(-6) 


-2.560(-6) 


0.95 


0.1635 


-9.213(-3) 


-1.034(-2) 


-9.786(-3) 


-2.393(-4) 


1.017(-3) 


2.096(-4) 


-5.569(-5) 


0.90 


0.3180 


-1.781(-2) 


-2.245(-2) 


-2.021(-2) 


-9.335(-4) 


3.996(-3) 


8.393(-4) 


-1.655(-4) 


0.85 


0.4623 


-2.593(-2) 


-3.643(-2) 


-3.147(-2) 


-2.033(-3) 


8.639(-3) 


1.841 (-3) 


-1.926(-4) 


0.80 


0.5947 


-3.403(-2) 


-5.209(-2) 


-4.389(-2) 


-3.470(-3) 


1.425(-2) 


3.050(-3) 


1.025(-4) 


0.75 


0.7137 


-4.334(-2) 


-6.855(-2) 


-5.794(-2) 


-5.159(-3) 


1.929(-2) 


4.050(-3) 


1.133(-3) 


0.70 


0.8177 


-5.704(-2) 


-8.293(-2) 


-7.450(-2) 


-7.021 (-3) 


2.019(-2) 


3. 860 (-3) 


3.615(-3) 


0.65 


0.9051 


-8.443(-2) 


-8.605(-2) 


-9.551(-2) 


-9.071(-3) 


7.270(-3) 


-4.744(-5) 


8.805(-3) 


0.60 


0.9768 


-0.1620 


-4.032(-2) 


-0.1273 


-1.192(-2) 


-5.589(-2) 


-1.615(-2) 


1.885(-2) 


0.55 


1.060 


-0.6220 


0.4026 


-0.2216 


-2.342(-2) 


-0.4951 


-0.1118 


3.530(-2) 


t 0.5244 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


0.50 


0.9322 


0.9441 


-1.289 


-1.317(-2) 


1.656(-2) 


1.043 


0.1764 


9.268(-2) 


0.45 


0.9589 


0.4818 


-0.8558 


-0.1228 


8.354(-3) 


0.5938 


4.776(-2) 


0.1726 


0.40 


0.9099 


0.4682 


-0.9012 


-0.1789 


1.439(-2) 


0.5780 


-4.266(-2) 


0.3572 


0.35 


0.8081 


0.5694 


-1.073 


-0.2463 


2.877(-2) 


0.6705 


-0.2372 


0.8163 


0.30 


0.6232 


0.8617 


-1.440 


-0.3617 


6.274(-2) 


0.9574 


-0.8477 


2.320 


0.25 


-0.2187 


3.306 


-3.834 


-0.9421 


0.3245 


3.571 


-7.001 


18.19 


t 0.2374 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


0.20 


0.8031 


-0.6934 


-0.2521 


-0.2627 


-9.644(-2) 


-0.8229 


3.568 


-9.625 


0.15 


0.4650 


-8.768(-2) 


-1.024 


-0.6262 


-2.999(-2) 


-0.2098 


2.168 


-6.409 


0.10 


0.2583 


7.056(-2) 


-1.555 


-1.133 


-1.064(-2) 


-6.378(-2) 


1.877 


-6.251 


0.05 


0.1006 


0.1021 


-2.758 


-2.463 


-2.574(-3) 


-1.361(-2) 


1.885 


-7.055 



As shown in Table II, the velocity field of 1PN order has an X3 component which 
is written as 

^3 = d 3 (p PN = 2sxix 2 x 3 , (7-7) 

where we use Eq. (2-41). This feature has also been found in the study of the 
Dedekind ellipsoid. v3 In the 1PN (and relativistic) case, we believe that the velocity 
field cannot be restricted in the planes orthogonal to the figure rotation axis in 
general. 

In Table III, the angular velocity, energy, angular momentum, luminosity of grav- 
itational radiation and velocity field of 1PN order in the ellipsoidal approximation 
are presented. There appear no singularities which are associated with instabilities 
because we do not take into account deformation of the ellipsoidal figure. From 
Table III, it is found that the velocity field of 1PN order has an X3 component, as 
in the case of the exact calculation. This implies that it is necessary to include the 
X3 component of the velocity field even in the ellipsoidal approximation. 

In Fig. 1, we show the angular velocity of Newtonian and 1PN orders as a 
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Table III. The 1PN angular velocity, energy, angular momentum, luminosity of gravitational ra- 
diation and velocity field in the ellipsoidal approximation along the equilibrium sequence of 
the irrotational Riemann S-type ellipsoid. The subscript e denotes the case of the ellipsoidal 
approximation. 



a 2 /ai 






JpNe 


{dE/dt) PNe 




<7c 


f e 




0.99 


8.050(-2) 


-4.284(-2) 


2.954(-5) 


5.154(-6) 


3.336(-2) 


-1.921(-3) 


-1.907(-3) 


-1.914(-3) 


0.95 


7.826(-2) 


-4.235(-2) 


7.275(-4) 


1.343(-4) 


0.1633 


-9.959(-3) 


-9.595(-3) 


-9.787(-3) 


0.90 


7.735(-2) 


-4.073(-2) 


2.849(-3) 


5.676(-4) 


0.3172 


-2.088(-2) 


-1.938(-2) 


-2.021(-2) 


0.85 


7.844(-2) 


-3.780(-2) 


6.259(-3) 


1.354(-3) 


0.4608 


-3.287(-2) 


-2.948(-2) 


-3.149(-2) 


0.80 


8.138(-2) 


-3.336(-2) 


1.084(-2) 


2.562(-3) 


0.5929 


-4.601(-2) 


-4.012(-2) 


-4.388(-2) 


0.75 


8.592(-2) 


-2.713(-2) 


1.645(-2) 


4.276(-3) 


0.7120 


-6.030(-2) 


-5.165(-2) 


-5.774(-2) 


0.70 


9.168(-2) 


-1.884(-2) 


2.296(-2) 


6.596(-3) 


0.8169 


-7.563(-2) 


-6.466(-2) 


-7.353(-2) 


0.65 


9.812(-2) 


-8.131(-3) 


3.027(-2) 


9.633(-3) 


0.9056 


-9.168(-2) 


-8.001(-2) 


-9.187(-2) 


0.60 


0.1045 


5.362(-3) 


3.835(-2) 


1.350(-2) 


0.9762 


-0.1079 


-9.899(-2) 


-0.1136 


0.55 


0.1099 


2.202(-2) 


4.729(-2) 


1.826(-2) 


1.026 


-0.1232 


-0.1234 


-0.1398 


0.50 


0.1130 


4.218(-2) 


5.749(-2) 


2.388(-2) 


1.054 


-0.1360 


-0.1560 


-0.1720 


0.45 


0.1126 


6.598(-2) 


6.993(-2) 


3.009(-2) 


1.055 


-0.1439 


-0.2003 


-0.2122 


0.40 


0.1071 


9.315(-2) 


8.678(-2) 


3.621(-2) 


1.026 


-0.1438 


-0.2611 


-0.2633 


0.35 


9.567(-2) 


0.1226 


0.1130 


4.096(-2) 


0.9658 


-0.1320 


-0.3441 


-0.3289 


0.30 


7.828(-2) 


0.1519 


0.1599 


4.251(-2) 


0.8710 


-0.1049 


-0.4562 


-0.4147 


0.25 


5.664(-2) 


0.1763 


0.2559 


3.900(-2) 


0.7430 


-6.093(-2) 


-0.6049 


-0.5289 


0.20 


3.434(-2) 


0.1886 


0.4789 


2.994(-2) 


0.5874 


-3.273(-3) 


-0.8009 


-0.6868 


0.15 


1.596(-2) 


0.1799 


1.078 


1.763(-2) 


0.4161 


5.746(-2) 


-1.070 


-0.9242 


0.10 


4.719(-3) 


0.1427 


3.088 


6.679(-3) 


0.2469 


0.1024 


-1.512 


-1.357 


0.05 


4.973(-4) 


7.597(-2) 


14.52 


9.189(-4) 


9.926(-2) 


0.1057 


-2.718 


-2.594 



function of the axial ratio aija\. Figures 1(a) and (b) are drawn for the Newtonian 
solutions, and (c) and (d) are for the 1PN solutions. The solid and dotted lines 
denote the results of the exact calculation and those of the ellipsoidal approximation, 
respectively. In the Newtonian case (or in the spherical case), the above two lines 
coincide. On the other hand, in the 1PN (and non-spherical) cases, the differences 
between the two lines appear, but are not large. 

The energy and angular momentum of Newtonian and 1PN orders are shown in 
Fig. 2 as functions of the axial ratio 02/01- In Fig. 2, we present the results only 
for 1 > 02/01 > 0.55, i.e., for the stable branch. The solid and dotted lines denote 
the same quantities as those in Fig. 1. The luminosity of gravitational radiation is 
shown in Fig. 3. The solid and dotted lines are again the same as in Fig. 1. From the 
figures, it is found that with increasing the 1PN correction, the angular velocity and 
angular momentum increase, while the total energy and luminosity of gravitational 
waves decrease. These features are reasonable, because with increasing relativistic 
correction, the self-gravity of the star becomes strong (and hence the binding energy 
becomes large), and to support the self-gravity, angular velocity must be larger. 

§8. Summary and discussion 

We have calculated an equilibrium sequence and the luminosity of gravitational 
waves of an irrotational and incompressible star in the 1PN approximation. We have 
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Fig. 1. The angular velocity of Newtonian (a) and (b) and 1PN (c) and (d) orders as functions of 
the axial ratio 02 /01. Figures 1(a) and (c) are magnifications of Figs. 1(b) and (d). The solid 
and dotted lines denote the results of the exact calculation and of the ellipsoidal approximation, 
respectively. 



developed a scheme to solve the hydrostatic equations for an irrotational fluid in the 
1PN approximation. By extending the present study, we will investigate irrotational 
binary stars in the 1PN approximation. Also, the 1PN solutions presented in this 
paper will be useful when examining the accuracy of numerical code for obtaining 
relativistic irrotational stars. 

In addition to establishing the method, we have an interesting result. The 1PN 
velocity field has an X3 component. This fact implies that internal motion never 
remains in the planes orthogonal to the figure rotation axis. (Note that such a waving 
motion closes even if it has an £3 component.) We can find in the Chandrasekhar and 
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1 0.9 0.8 0.7 0.6 1 0.9 0.8 0.7 0.6 



(a) (b) 




1 0.9 0.8 0.7 0.6 1 0.9 0.8 0.7 0.6 

a 2 / a, a 2 / a ( 

(c) (d) 



Fig. 2. The total energy of Newtonian (a) and 1PN (c) orders, and the total angular momentum of 
Newtonian (b) and 1PN (d) orders as a function of the axial ratio a 2 /ai. The solid and dotted 
lines are the same as in Fig. 1. 

Elbert paper© that there is also an £3 component of the internal motion in the case 
of the 1PN Dedekind ellipsoid. Moreover, Uryu and EriguchiElP numerically showed 
that the internal flow of the compressible Dedekind-like star has an X3 component 
even in the Newtonian case. Therefore, except for incompressible, Newtonian cases, 
non-axisymmetric stars with an internal motion will have an X3 component of the 
velocity field in general. 

Finally, we discuss the validity of the ellipsoidal approximation. One can see 
from Figs. 1 ~ 3 and Tables I ~ III that the error of the 1PN angular velocity, 
energy, and angular momentum in the ellipsoidal approximation are less than about 
~ 5% in the range of 1 > 0,2/0,1 > 0.7. Although some of the coefficients of the 
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Fig. 3. The luminosity of gravitational radiation of the Newtonian (a) and 1PN (b) orders as 
functions of the axial ratio aijax. The solid and dotted lines are the same as in Fig. 1. 



1PN velocity potential, q and f, deviate from exact values by a large factor near 
(12/ai ~ 0.55, the error is not so large for larger 02/01. Therefore, we may conclude 
that the ellipsoidal approximation is a good one if the star is not highly nonspherical 
(1 > a%ja\ > 0.7). This result is encouraging for study of binary stars in the 
ellipsoidal approximation becausein the binary case, it is expected that 02/01 and 
03/01 do not become very small. 
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